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Abstract 

A Gibbs-like approach for simultaneous multi-scale correlation 
functions in random, time-dependent, multiplicative processes for 
the turbulent energy cascade is investigated. We study the optimal 
log-normal Gibbs-like distribution able to describe the subtle effects 
induced by non-trivial time dependency on both single-scale (struc- 
ture functions) and multi-scale correlation functions. We provide an- 
alytical expression for the general multi-scale correlation functions 
in terms of the two-point correlations between multipliers and we 
show that the log-normal distribution is already accurate enough to 
reproduce quantitatively many of the observed behavior. The main 
result is that non-trivial time effects renormalize the Gibbs-like ef- 
fective potential necessary to describe single-time statistics. We also 
present a generalization of this approach to more general, non log- 
normal, potential. In the latter case one obtains a formal expansion 
of both structure functions and multi-scale correlations in terms of 
cumulants of all orders. 

Key words: Turbulence, Multifractals, Stochastic Processes, Shell 
Models 

1 Introduction 

Small scales, three-dimensional, turbulent fluctuations are sustained by 
the energy cascade mechanism: energy is injected at large scales, Lq, and 
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dissipated at small scales, r\. The statistical properties of the energy trans- 
fer throughout the range of scales going from Lq to r\ (inertial range) are 
thought to possess highly non-trivial features, mainly, but not only, summa- 
rized by the presence of anomalous scaling, i.e. pth order velocity structure 
functions, S p (r) = ((v(r) — v(0)) p ) have a power law scaling with anomalous 
exponents, S p (r) ~ (r/£o)^ ■ The existence of a net, direct, energy flux 
in the inertial range is a clear indication of the out-of-equilibrium nature 
of turbulent flows. The exact symmetry of the inertial Navier-Stokes terms 
with respect to direct or inverse energy transfer is explicitly broken by the 
existence of an energy source at large scales and an energy sink at small 
scales. The energy cascade process has been often, and fruitfully, described 
in terms of a multi-step fragmentation process describing the tendencies of 
inertial range eddies to break in smaller and smaller eddies, following the 
celebrated Richardson scenario Q . The spatio-temporal complexity of the 
fragmentation process has been successfully described by Q using the mul- 
tifractal language, which have proved able to reproduce qualitatively and 
quantitatively single-scale, multi-scale and multi-time multi-scale velocity 
correlation functions ||, || || . 

By using the multifractal language one can assume that the velocity differ- 
ence on scale r 2 is linked to the velocity difference at scale r\ > r 2 by the 
equation 

Sv{r 2 )=M{r 2 ,r 1 )5v(r 1 ), (1) 

where M(r 2 ,ri) is a suitable random variable. Moreover, one often limits 
the complexity by assuming almost uncorrelated multipliers: 

M(n,r 3 ) = M(n,r 2 )M(r 2 ,r 3 ) 

for any choices of r\ < r 2 < r 3 in the inertial range. The above assumptions 
lead to the concept of random multiplicative process as a possible way to 
characterize the turbulent, multifractal, velocity field. In ||] a Gibbs mea- 
sure for the random multiplicative process have been introduced in order 
to compute, self consistently, the anomalous scaling exponents for a class 
of Shell models. Recently, the idea discussed in (|] have been further de- 
veloped in by studying some new proposals on how to use a Gibbs-like 
approach to describe energy fluctuations in the inertial range. According 
to the Kolmogorov hypothesis of local interactions in Fourier space, one ex- 
pects short range correlations between multipliers, while velocity fields may 
still show long-range correlations. On one hand, the simplest phenomeno- 
logical multifractal description able to capture the correct anomalous scal- 
ing for single-scale structure functions is based on the assumption of a 
complete uncorrelated multiplicative process ||. On the other hand, the 
assumption of completely independence between multipliers in untenable 
because in disagreement with some theoretical and experimental results. 
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In this paper, we try to perform a step toward more realistic description 
of turbulent energy cascade by studying also deterministic and stochastic 
processes with correlation between multipliers. Correlations are introduced 
by a non-trivial time evolution of multipliers in a self-consistent way, i.e. 
by imposing that time and spatial fluctuations are linked as dictated by 
the non-linear terms of the Navier-Stokes eqs. The main goal is to ob- 
tain the optimal effective Gibbs potential able to describe the simultaneous 
fluctuations of the velocity field at all scales. 

We study the problem in both deterministic models of turbulent energy 
cascade (Shell Models) fi |o) and in stochastic, time-dependent, multi- 
plicative processes. Let us recall that shell models are the simplest deter- 
ministic models whit anomalous multi-time and multi-scale velocity corre- 
lation functions. Moreover, shell models posses the special feature to have 
no sweeping terms, i.e. also temporal properties of the energy cascade can 
be tested. 

The paper is organized as follows. First we briefly recall the main 
features of Shell models and we define the set of observable we want to 
describe within the Gibbs approach. Then, we discuss the simplest, log- 
normal, approximation for the Gibbs-potential. Log-normal potential en- 
joys the properties to be "exactly solvable" , any structure functions and 
multi-scale correlation function possess an explicit expression in terms of 
the two-point spin-spin correlation function (see below). We present some 
evidences that already the simple log-normal approximation is able to cap- 
ture many of the observed behaviors for both stochastic process and shell 
models. Further, we discuss how to generalize the approach to a more gen- 
eral, i.e. non log-normal, potential, and finally, we conclude with comments 
on possible future works. 



2 The Gibbs ensemble 

Shell models describe the energy turbulent transfer on a set of scales (shells) 
in the Fourier space, k n = ko\ n , where fco is the smallest wavenumber and 
A is the inter-shell ratio, usually set to 2. Velocity shell variables, u n (t) 
axe complex numbers representing velocity fluctuations, S rn v, over a scale 
r n = fc" 1 , with n = 0,...,N. Among all possible shell models a very 
popular one is the Sabra model jfl| an improved version of the GOY model 
[§. The model is: 

(d/dt + vkl)u n = ik n (u n+2 u* l+1 + bu n+ iu* l _ 1 - cu n _iu„_ 2 ) + /„ (2) 

where v is the kinematical viscosity, 6, c are free parameters fixed by requir- 
ing that energy and helicity are inviscid invariants of the models and /„ is 
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the external forcing with supports only on large scales (small shell indexes). 
The existence of only local interactions between shells allows to have highly 
non-trivial time properties, i.e. the model is a reliable approximation of 
velocity evolution in a quasi-Lagrangian reference frame. The model is 
known to possess realistic multi-time and multi-scale correlation functions, 
including anomalous inertial range scaling and dissipative anomaly, for a 
review see || |l |l(|. Typical observable checking single-scale Probability 
Density Functions (PDF) are given by structure functions: 

S P (k n ) = (\u n \P). (3) 

On the other hand, multi-scale single-time observable which will play a 
relevant role in the following are the two-scale correlation function: 

F p<g (k n ,k n+m ) = (\u n \ p \u n+m \ q ). (4) 

Usually, multifractal phenomenology is based on uncorrelated multipli- 
ers, M{k n ,k n+m ), connecting the two shell velocity at scales k n ,k n + m as 
Un+m = M{k n ,k n+m )u n . This simple scenario leads naturally to pure, 
anomalous, power law scaling for the structure functions (||), 

Sp{kn) ~ k n ^^ 

and to pure fusion-rules predictions for the two scales correlation [Q, ^, |E| : 

/fc„)-«5)fc-C(f+?). (5) 

The numerical integration of shell model equations (^) shows that the be- 
havior of the multi-scale correlation function, F p . q {k n , k n +m)i predicted by 
(||) is true only asymptotically for very large scale separation k n + m /k n — * 
co, while important deviations are detected for scale separation, k n+m /k n ~ 
0(1). The origin of such deviation can easily be understood if one looks 
at the typical temporal evolution of the energy contents at different scale. 
In shell models, energy is transferred down-scale by a burst-like activity, 
strong coherent energy bumps travel from large-scales to small-scales. Each 
scale has its typical, fluctuating, eddy-turn-over time, t„ ~ l/(u n k n ). En- 
ergy is transferred from scale k n to scale k^ in a typical time r„ jv = t„ — tjv 
(with N > n). The non-instantaneous, intermittent, propagation of energy 
from shell to shell is dynamically realized by non-trivial fluctuations of shell 
model phase variables, <f> n of (||) - where u n — \u n \ exp 1 ^™. 
The time-delay in the information propagation has important feedback also 
on single time observable as structure functions and multi-scale correlation 
functions (see below) . This is the physical reason why the understanding of 
single time statistics calls for the understanding also of multi-time statis- 
tics. The scope of this article is to investigate to which extent one may 
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try to build up an effective Gibbs-like description capable to incorporate 
the effects of non-trivial time fluctuations on single-time statistics. In the 
following we restrict ourself to discuss the probability distribution function 
of shell amplitude, |u n |. In jffl, following the original proposal made in 
a Gibbs hypothesis has been developed for the simultaneous probability 
distribution function of the u n set of shell variable: 

P( Ul , M2 ,..., Ujv )oce-*^<" 2 --<"~> (6) 

where for sufficiently high Reynolds numbers we will suppose the potential, 
to become translational invariant in the shell index (independent on 
the UV and IR boundary condition). As previously said, we also make, 
following j|, ||] , the further assumption that the potential depends only on 
ratios between shell variables, i.e. it is an homogeneous function of zero 
degree. The physical rationale for this assumption stem from the original 
Kolmogorov remark that energy-cascade is maintained by a local transfer 
in Fourier space. Introducing the spin variables, <Xj = log 2 (|ttj l/lifj+i |) we 
may rewrite the Gibbs-hypothesis for the shell amplitudes as: 

$(|«i|, \u 2 \, ■ ■ • , |«jv|) = 0(oi,er 2 , ■ ■ • ,0- A^- 
Let us remark that the Gibbs-potential has nothing to do with the "equilib- 
rium" distribution obtained by the equipartition of the inviscid invariants 
even in the limit of v — ► 0. In R| a very detailed numerical and theo- 
retical analysis on the consequences of such a description has been made, 
giving strong evidences that the formalism is consistent with the statistical 
properties of the model. Here, we want to push further this hypothe- 
sis by explicitly looking at the possible Gibbs-potential able to reproduce 
quantitatively and qualitatively the measured structure functions (^) and 
multi-scale correlation functions (^). A pure uncorrelated multiplicative 
process on the shell amplitude is described in the Gibbs formalism by a 
simply infinite temperature Gibbs potential with no interaction between 
spins: 

JV 

tf>(ai,a 2) ... ) (TN) =^2V((Tj). (7) 
i=i 

Indeed, it is simple to show that in this case we have for the structure 
functions the exact anomalous scaling: 

/JV n 
J] da, exp{-(]T V{aj) +plog(2) £ a,)} « fc-«*> 

with = ^o<72(2 _p<T ), where with ( • } we mean the average with respect 
to the Gibbs measure Y[ d&j exp{— 4>(cr±, 02, . . . , Ojv)}. Similarly, for the 
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multi-scale correlation function we have: 

F p , q (k n ,k n+m ) = (\u n \P\u n+m \ q ) = <2~ (p+9) ^i CTj ~ ? £"-+^), (8) 

which correspond to the fusion-rule prediction (|^) for all scale separation. 
On the other hand, as previously said, the pure uncorrelated hypothesis 
cannot be correct because it does not predict the observed behavior of 
multi-scale separations for small scale separation, k n /k n+m ~ 0(1). In- 
deed, two-scales correlation as (^) tends to the fusion-rule prediction (||) 
only asymptotically -see for example fig. (|^). Thus, it seems necessary to 
add to the uncorrelated potential (0) some spin-interaction. The simplest 
analytical way to do it is to stay within all possible interaction in a Gaus- 
sian field, i.e. to consider a "correlated" log-normal distribution for the 
shell variables. As we shall see in the last section, this assumption is not 
restrictive, most of the qualitative and quantitative results here presented 
can be extended to more complex probability distribution. 

A log-normal uncorrelated process is described by the Gibbs potential: 

02, . . • , <Jn) — J2f=i (2c?T wnere the only two free parameters are 
the log-mean ho and the log- variance c 2 . Correlations among multipliers 
can be introduced by writing the Gibbs potential as: 

</>(<7l,<7 2 , . . . ,<7jv) = T^'^VjAjiCTi - h /c 2 <7j (9) 

3,i 3 

where now the matrix of interaction is given by Aji — Sji / c 2 — Jji . Clearly, 
by taking Jji depending only on the separation between shell indexes, j — i, 
we may describe the most general translational invariant, log-normal, two- 
body potential. Translational invariance in the spins variable is the coun- 
terpart of scaling invariance for the velocity variables. 
As previously said, in shell models one cannot simply disentangle the ampli- 
tude fluctuations, from the phase fluctuations, (f> n , in other words one 
cannot expect to reproduce quantitatively the multi-scale fluctuations with- 
out a further explicit introduction, in the Gibbs formalism, also of phase- 
variables. In order to make the discussion simpler, we will test quantita- 
tively the Gibbs-formalism on a stochastic, time-dependent, multiplicative 
signal involving only amplitude fluctuations, meant to mimic the amplitude 
evolution of shell models dynamics. We will go back to comparison with 
the outcomes of the original deterministic dynamics @ only to show the 
ability of the Gibbs formalism to catch the main qualitative behaviors. 
The stochastic, time-dependent, multiplicative process is built as follows 
(see also |||). 

We introduce N i.i.d. random variables, Wj = M(fc J+ i, kj), one for each 
shell, describing the uncorrelated instantaneous multipliers connecting am- 
plitudes of shell variables between shells n and n + 1, i.e. |ii n +i| = W n |u n |. 
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The probability of W coincides with the log-normal uncorrelated Gibbs- 
potential: P(W n ) oc exp{— 2^ }• ^° g enera te the time dynamics we 
proceed as follows. We extract W n with probability P(W n ) and keep it 
constant for a time interval + r„], with r n = l/(|u n |fc„) being the local 
instantaneous eddy-turn-over time, Thus, for each scale k n , we introduce 
a time dependent random process W n (t) which is piece- wise constant for 
a random time intervals \t^n \ t^ + r„], if is the time of the kth jump 
at scale n. The corresponding velocity field at scale n, in the time interval 

(k) (k) 

t n < t < tn + T n , is given by the simple multiplicative rule: 

\u n (t)\=W n (t)\u n ^(t^)\. (10) 

What is important to notice is that at each jumping time, t„ , tn , ■■, t^\., 
for any scale, n, only the local velocity field is updated, i.e. information 
across different scales propagates with a finite speed. It is easy to realize 
that the propagation speed is proportional to the characteristic speed of 
the energy cascade in turbulent flows, i.e. a fluctuations in the multiplier 
at scale k n takes a time T ~ r„ — t„, to propagate down to scale k m . In 
this way we reproduce the phenomenology of the non-linear evolution of 
the shell model dynamics: d/dt „, which is itself meant to mimic 

the non-linear evolution of a Navier-Stokes field in quasi-Lagrangian ref- 
erence frame. This is the simplest stochastic evolution with non-trivial 
spatial and temporal fluctuations in qualitative agreement with the shell 
models phenomenology In the following, whenever we refer to the 

time evolution of the stochastic process at scale k n we use the notation 
Un^ (t) to distinguish it from the time evolution of the deterministic shell- 
model velocity, u n (t). The link between spatial fluctuation and temporal 
fluctuations, r„ = l/(k n u n ), has important feedback even on the single 
time structure functions, S p (k n ), and single time multi-scale correlation 
functions, F Ptq (k n ,k n+m ). In fig. |l| we plot the scaling exponents of the 
structure functions measured on the stochastic field, S p s \k n ) = (\un^\ p ) 
with and without fluctuating eddy-turn-over time. As one can see there is 
an important "renormalization" effect when eddy-turn-over times fluctu- 
ate. The scaling exponents move from the simple uncorrelated value when 
the time dynamics is trivial to a renormalized value when local multipliers 
are updated with stochastic times. This renormalization effect can be un- 
derstood in terms of the correlation between the fluctuating eddy-turn-over 
times and the multipliers |l^] . This is the first evidence of a relevant effect 
on the energy cascade introduced by the time dynamics even on single-time 
observable. In fig. || we also compare the behavior of a single time multi- 
scale correlation function, F2^(k n , k n+m ), measured on the stochastic pro- 
cess with and without time dependency and in the original shell model (^) . 
It is important to notice that both shell model and the stochastic, time- 
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dependent, process show the same similar slow approach to the asymptotic 
plateau for the normalized multi-scale correlation function. The above re- 
sult shows that the departure measured for small scale separation from the 
asymptotic fusion-rule prediction (^) is mainly due to non-trivial correla- 
tion between multipliers introduced by the time-dynamics 0. 
In fig. [3] we show the spin-spin correlation, i.e. the correlation among mul- 
tipliers, for both the time dependent random multiplicative process and 
the equivalent quantities in the shell model. It is rather clear that time 
dynamics introduces a correlation among multipliers in a non trivial way. 
Moreover, it is important to remark that the spin-spin correlation is qual- 
itative similar for both the time dependent random multiplicative process 
and the shell model. In particular let us remark that the near neighborhood 
correlation is negative in both cases. 



3 Some rigorous results 

In both shell models and the stochastic process all single time correlation 
functions arc determined by the whole spatio-temporal dynamics. We want 
now investigate the possibility to reproduce the effects of the non-trivial 
time properties of the energy cascade on the single-time statistics . We try 
to do it by an effective, time independent, Gibbs potential. Therefore, the 
Gibbs-potential must be seen as a renormalized set of time-independent 
interaction describing the whole set of possible multi-scale single-time cor- 
relation functions. 

In order to keep our discussion as simple as possible we confine here 
to investigate log-normal distributions. In the next section, we generalize 
our results for any probability distribution of random multipliers. As dis- 
cussed in section 2, we define the random multipliers as the ratio between 
stochastic fields at neighboring scales: u^^/u^ — Wi = 2~ ai . It follows 
that the joint probability distribution for a given set of spins/multipliers 
{a.j} is given by: 

N 1 

P[ai] oc exp{/i /c 2 ^2 a ' 1 ~ 2 51 1 ' -' rr ' rr ' ' • ( n ) 

i=l i,j 

where Aij is the spin-spin interaction. We introduce the shorthand nota- 
tion a to denote the A-component vector formed by all spins a — (a\ , 02, • • • , 
and A to denote the interaction matrix. Then, we may rewrite the partition 
function Z(ho): 

Z(h ) = I daexp{-^aAa + h Q ■ a} ~ exp^ oi ~ 1Ko) (12) 
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where the vector ho is made of constant entries for all scales: ho = ho/c 2 (1 , 1 , 1 , 1 , 
Having restricted ourselves to the most general Gaussian distribution it is 
not surprising that one can work out an explicit formula for the most gen- 
eral multi-scale correlation functions in terms of the two-point connected 
correlation function, {<Ji(Jj) c = (ci<7j) — (<Ti)(crj), only. In the appendix we 
present the long straightforward calculation leading to the final expression: 



r ,i,„! ^ \"--»»»+m)«) = r T;; 



"P,1 - ■',// • .„/-</•, ~ e P,q ( 13 ) 
with 

n 

T™ q = (log2) 2 pq(^2{<JiCr n+ i) c + (ai<7 n+2 ) c H V (<Ti<T n+m ) c ). (14) 

i=l 

The above formula has a particularly appealing interpretation: deviations 
of the multi-scale correlation functions from its "multiplicative uncorre- 
cted" fusion-rules prediction (||) is governed by the short range correlation 
between multipliers. Indeed, in the RHS of Jl^) the main contribution is 
carried by the first connected correlation function, (a n a n ^i) c , while all the 
other terms becomes less and less important because they connect spins at 
larger and larger distances. It turns out that spins are anti-correlated at 
short distances (see fig. ^). This is the reason why all normalized multi- 
scale correlation functions (^) converge to a plateau smaller than unity, 
i.e. Cp% < 1. Moreover, in this log- normal approximation, the coefficients 
Cp™^ are symmetrical in p, q something which may be exploited in order to 
reduce the number of degrees of freedom in closures jl3| . 



3.1 Numerical Tests 

Let us now try to check numerically whether the log-normal approximation 
is in good agreement with the numerics observed in the stochastic signal, 
Un\t). As already remarked, the most general log- normal distribution is 
completely fixed once one defines the coupling matrix and the magnetic 
field ho defining the "bare" probability distribution function of multipliers, 
ho = — (log 2 W). Moreover, the coupling matrix is in one-to-one corre- 
spondence with the connected two-spins correlation functions: 

*V = Mo- (15) 

We have therefore taken as our best guess for the coupling potential the 
expression for Jy obtained from (|l5|) by using the measured two-point cor- 
relations (cjOj) in the stochastic evolutions -a plot of (ciCTj) c is shown in 
fig. ra. Than, from this numerical input we can check whether the quantities 
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analytically computable from (O) are in agreement with those measured on 
the stochastic signal. In fig. W we compare the two multi-scale correlation 
functions calculated either numerically or from the ( |ll| ) expression. As on 
can see the agreement is qualitative and quantitatively very satisfactory. 
The above results tell us that the most general log-normal distribution, 
chosen to exactly reproduce the measured two-point correlation function, 
(cTjOj) c is also able to reproduce in a quantitative way the multi-scale corre- 
lation function with high accuracy, i.e. the log-normal approximation with 
the potential given in fig. ^ is a very close approximation of the "effec- 
tive" single-time probability distribution of the complete time-dependent 
stochastic process. Also scaling exponents measured on the correlated log- 
normal potential are in good agreement with those measured numerically 
on the time-dependent stochastic process (with deviations of order 5% on 
the 10-th order exponent). Furthermore, in the next section, we will show 
that the log-normal result can be seen as the first term in a systematic 
expansion in cumulants of the most general probability distributions. 



4 A generalization to non-Gaussian distribu- 
tions 

We want here to present a simple argument showing that the expression 
( p"4| ) obtained within the log-normal approximation can always be seen as 
the first term of a formal cumulant expansion for the most general poten- 
tial. 

Let us consider the generic interacting potential $(cti, ...,ctjv) among the 
N spins. Where we now may include also three-body and multi-body in- 
teractions in $ and/or the spins variable can take values also on a discrete 
set (Ising-like systems) . If we go back to the observable we want to control 
in our multiplicative process, we realize that one may write both struc- 
ture functions and multi-scale correlation functions as suitable partition 
function calculated with suitable external magnetic field: 

«) oc ex P{ ? ' Hn - Qipu ctn)} = Z0D (16) 
{*<} 

and 

(«+m) « ex P{^ ' H^ q n+m - $((7!, ...,a N )} = Z(H™ +m ) (17) 

where and H^ L+m are site-dependent magnetic fields which are explic- 
itly defined by the expressions ( |26|) a nd ( pg)in appendix A -see appendix 
A also for a detailed derivation of (|16|) and (17]). Notice that in the generic 
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potential <& we already include any possible linear coupling with an homo- 
geneous magnetic field. Exploiting the expansion in cumulants of a generic 
partition functions one easily obtain the formal expression for the structure 
function and the two-scale correlation function, respectively: 

log«)oc^i (^^■■■^ k )cH? l (i 1 )HP( l2 )---Hr(i k y i (18) 

k i 1 ,i 2 ,---,ik 

^og{u p n u q n+m ) ( CT n CT i2---^Jc^;* +m (ji)ff^ + „ i («2)---ff^ +m (ife), 

k »x,ij 

(19) 

where all connected correlation functions, (t7 il cr i2 ...er iii } c are calculated at 
zero external magnetic fields, H — 0. It is easy now to realize that the 
general expression for, Cp™\ i.e. the deviation from the pure fusion-rules 
prediction for two-scale correlation function can be expressed as a power 
series of suitable combination of external magnetic fields: 

log(CW)=^i Y, (^^•••^>cM|:i,4 ( 20 ) 

k i 1 .i 2 ,---,ik 

with 

= II ( H Z+mU) + Hl{ 3 ) Hl +m { 3 ) - . (21) 

j=ii,...,i k 

Obviously, for a given order, fc, the different indexes ii,...,ik must take 
values between 1 and n + m; indeed, it is sufficient that only one among 
ii, ik does not fall between n and n + mto have that the external mag- 
netic fields HP, HI, H^ +m , H^ +m vanish and therefore Mf ,q i = 0. 



Expression (20) tell us that the previous log-normal result ( [14| ) can be seen 
as the first contribution, k — 2, in the above expansion, being the contri- 
bution for k = 1 identically zero. It is easy to realize that if the potential is 
exactly log-normal, only contribution form the k = 2 term appears, while 
in the most general case one need to control also three-point correlations, 
(ercer ) and multi-point correlations. 

In appendix B we develop explicitly the expression of ( pcf ) for any order k 
and also the similar interesting expansion obtainable for single scale struc- 
ture functions dig). It is important to remark, that from the latter one 



also obtain a formal expansion of £(p) exponents in power of p: 

c(p) = 5>y (22) 

3>0 

where the set of Cj are connected to the choice of the interacting potential, 
$. In the case of a log- normal potential we nay write the explicit expression 
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of scaling exponents in the limit of small enough scales, i.e. for n large 
enough: 

c(p) = p- £ fo>c - £«<™> c + ( CT ^i> c • • • (^i) c )- (23) 

4=1 4=1 

The above expression tells us that the renormalization of exponents due to 
the appearance of correlation among multipliers is linked to the magnetic 
suscettivity, i.e. to the integral of the connected two-point correlation for 
all scale separations. Obviously, in the simple case of pure uncorrelated 
log-normal process with (uiUj) c — Si.j ■ c 2 and (<Ji) c = ho, one finds the 
well known result: 

C(p)=p/i -ip 2 c 2 log(2). (24) 

5 Conclusions 

A Gibbs-like approach for single time multi-scale correlation functions in a 
class of random multiplicative process with non-trivial time dependencies 
has been investigated. We have shown that there exists an optimal log- 
normal Gibbs-like measure able to reproduce with high accuracy the effects 
induced by the temporal dynamics on the single-time correlation functions. 
We have explicitly calculated the expression of the most general two-scale 
correlation functions (u^u^ , m ) in the log-normal approximation. We have 
also shown that the log-normal result can be seen as the first order of 
a formal cumulant expansion obtained for a completely general potential 
of interactions between spins (multipliers). Within this formalism, also 
scaling exponents, £(p) have a simple power-law expansion in terms of the 
order of the moment, p. The first two terms in this expansion coincides 
with the usual quadratic log-normal expression. 

Qualitatively, we expect that very similar results can be obtained in order 
to describe the multipliers statistics of more realistic models as the case of 
shell models. The qualitative similar behavior shown in figs. H|| indeed is a 
good evidence that the stochastic process here studied mimics quite well the 
shell model dynamics. The presence of non-trivial phase-phase correlations 
and phase-amplitude correlations in shell models is the major difficulties to 
overcome if one wants to apply in a quantitative way the Gibbs approach 
in this case. To proceed on this route one needs to introduce some spin 
variables describing phase fluctuations and some new potential describing 
phase-phase interactions and phase-moduli interactions (see also 0). Also, 
in shell models higher than second order connected correlation functions 
appears, as can be easily checked numerically, and therefore the log-normal 
approximation must be meant only as the first order term in the cumulant 
expansion as previously discussed. 
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5.1 Appendix A 

In order to compute analytically the expression for any correlation func- 
tion, it is useful to realize that the calculation of either structure func- 
tions or multi-scale correlation function can be reduced to the calcula- 
tion of a particular partition function with a suitable non-homogeneous 
magnetic field, Z(H ,p,q), where with H we intend the one-dimensional 
vector whose N components are given by the magnetic field in each site: 
H = (hi, h%, ■ ■ ■ , Jin). In particular, from ( p"l| ) and (||) one gets for the 
structure function: 

N 

JJ^exp^^i^-P^WE^^-iE^^.*^ . (25) 

i=l 

The extra term in the exponential plog(2) Yn=i ai can ^ e seen as an a< ^" 
ditional, position dependent, magnetic field of intensity p. Therefore, the 
structure function (u p n ) is proportional to the original partition function 
with a modified magnetic field 

«} oc Z(Hl) oc exp^S+^^S+H 

where is a vector with components given by: 

fl»(i) = -0(n-i)plog(2), (26) 

where we have introduced the Heaviside function, 9(x). In this notation, 
the non-homogeneous magnetic field acts from the integral scale i = 1 up to 
the inertial scale i = n. It is easy to realize, that similar expressions can be 
derived for the most general multi-scale correlation functions. In particular, 
the two scale correlation function, (u^it^ +m ) can be also calculated via a 
new partition function with a modified, site-dependent, magnetic field: 

«< +m ) cx Z{H^ n+m ) ~ exp^:^+^-^r,+ m +M, (27) 

where with H v n q n + m we denote the magnetic field vector whose i-th com- 
ponent is given by: 

fl£JUm(») = -(*(« - i)V + 0(n + m- i)q) log(2). (28) 

Now, it is long but simple to show with algebraic manipulation of previous 
expressions, that the prefactor, Cp$ defining the deviation of {u p n u!; l+m ) 
from the exact fusion rules prediction is given by: 
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with 

n 

T™ q = (log2) 2 pqC^2{<7 l cr n+1 } c + (o- J cr„ +2 }c H 1- (<Ti<r n +m)c), (30) 

i=l 

as reported in the text. 
5.2 Appendix B 

Let us here analyzed in more details the expression given in the body of 
the paper for the normalized two-scale correlation function: 

log(^=E^ E K^"-^>cMM ; ... ifc (31) 

with 

= II ( H Z + m(j) + Hl{j) Hl +m { 3 ) H?+i(j)) . (32) 

j=h,...,i k 

It is easy check that for k=l we may distinguish two cases (i) if 1 < i\ < n 

M™ = log2(-(p + q)-q + q+(p + q))=0; (33) 

(ii) if n < i\ < n + m: 

Mf' q = log2(-q- + q + 0) = 0. (34) 

Therefore we may conclude that the first order contribution in the cumulant 
expansion is identically zero. 

For k=2, the only non-vanishing contributions are those with 1 < ii < n 
and n < «2 < n + m, which give: 

Mf;% = (log2f{q{p + q)- q 2 ) = (logifpq (35) 

and after a permutation between %i and z 2 we have: 

n m 

log(C^)| fe=2 = (/o. 9 2) 2 TO ]T]>> 4( 7„ + ,} c ( 36 ) 

i=l 3=1 

which coincide with the log-normal contribution. 

For orders larger than 2, things becomes more complex, for example for 
k=3 we have two cases with non-vanishing contributions: 

• 1 < ii < n and n < ii, is < n + m leading to: 

M?f iaii3 = (log2f{-(p + q)q 2 + q 3 ) = -(log2) 3 (pq 2 ) (37) 
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• 1 < H>i2 < n an d n < < n + m leading to: 

■\ I = {log2f{-q{p + qf + q 3 ) = {log2)\-qp 2 - 2pq 2 ). (38) 

By considering all indexes permutation we get: 



n rn m 



log(C^)| fc=3 = -(W2) 3 (w 2 EEE^ (7 »+^«+^ 

i— 1 j—1 s—1 



n n m 



(2pq 2 + qp 2 )Y,J2Y,^ a i (7n + s ^y ( 39 ) 

2—1 j — 1 S — 1 

It is now possible to explicitly write down the generalization to any order 
k. If we consider the indexes ii, ...ik there are n\ indexes between 1 and n; 
ri2 indexes between n and n + m; with obviously n\ +ri2 = k. We therefore 
have: 

i°g(c&s>)u= E (-^2)*{[(p+gr(+«r-(9)*]- 

rii+n 2 — fe 
n n m rn 

E •- E E •- E ( c7 Ji- a 'in 1 ' S i+n-0- S „ 2 +n}c} (40) 

with ni > and 712 > 0. 

Let us now discuss how to obtain the expansion of scaling exponents 
in power of p. Starting from their definition: 

«> kn = T (4i) 

we have, using (|l9|): 

cw = -^Em E K--Oc^,... A (42) 

where: 

G l,.,i h = II ^(i) = (-^2) ft , (43) 

j=ii! — ,»* 

because all indexes ii, must be between 1 and n otherwise ^ = 

0. Therefore, we obtain the expression cited in the body of the paper: 

C(p)=J2c 3 P 3 (44) 
i>o 

with Cj defined in terms of the microscopic potential $. 
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Figure 1: Comparison between the scaling exponents, ((p), calculated on 
the simple multiplicative process without time dependencies, (x) and with 
fluctuating time, +. Notice the renormalization observed in the values of 
the exponents once fluctuating eddy-turn-over times are switched on [12] 
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Figure 2: Normalized multi-scale correlation function C, 9 = , " "+7 »+I 

> ("„+„,) ( M S ) 

with n = 12 as a function of the scale separation to, for shell model (+) and 
stochastic multiplicative signal with time dependencies (x). The straight 
line of value 1, corresponds to the trivial case of a multiplicative uncorre- 
cted process without any time-dependency. 
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Figure 3: Two point normalized connected correlation function among mul- 
tipliers, spins, {<Ji<Jj) J {(JiOi) c with i = 8 at changing the scale separation 
j = i — 8, . . . ,i + 8 calculated on the shell model (x) and on the time- 
dependent stochastic signal (+). Both i and j are in the inertial range for 
the shell model simulation. 
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Figure 4: Comparison between two-scale normalized correlation functions, 
calculated from the log-normal Gibbs formalism and on the stochastic time- 
dependent multiplicative process. Top curves: Cffi from the Gibbs for- 



malism, (+), and for the stochastic process, (x). Bottom curves: for 
the Gibbs formalism, (squares), and for the stochastic process (*) 



